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ABSTRACT 

The  practical  and  rigorous  solutions  of  the  potential 
problem  associated  with  the  harmonic  oscillation  of  a  sub- 
merged rigid  body  of  arbitrary  shape  is  presented.   The  use 
of  Green's  function  reduces  the  determination  of  the  poten- 
tial to  the  solution  of  an  integral  equation.   The  integral 
equation  is  solved  numerically  and  the  dependency  of  the 
hydrodynamic  quantities  such  as  added  mass  and  damping 
coefficients  of  the  object  on  the  frequency  of  the  oscilla- 
tion is  established. 

Several  checks  are  made  on  the  numerical  results.   These 
include  the  Haskind's  relations  check,  an  energy  check  for 
the  radiation  problem,  and  comparisons  with  closed  form 
solutions.   All  the  checks  and  comparisons  are  successful  and 
it  appears  that  the  numerical  procedure  employed  yields  valid 
and  accurate  results. 
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NOMENCLATURE 

Symbol  Description 

a  characteristic  length  of  object 

a  wavelength  parameter,  a  =  Zira/L,  dimensionless 

C.  wave  force  or  moment  coefficient  in  the  i 

direction,  dimensionless 

ds  elementary  dimensional  surface  area 

ds  elementary  dimensionless  surface  area 

E.  energy  transmitted  over  one  period 

f.,f  distribution  function,  dimensionless 

¥.  wave  force  or  moment  component  in  the  i 

direction,  dimensionless 


1 


g  acceleration  of  gravity 

G  Green's  function 

h  water  depth 

h  relative  depth  of  water,  dimensionless 

h.  functions  defined  by  (19),  dimensionless 

H  wave  height 

Im  imaginary  part 

n  ,n  ,n  unit  normal  vectors,  dimensionless 

x*  y'    z                      ' 

P'  dynamic  pressure  associated  with  the  wave 
interaction  with  the  fixed  object 

P.  dynamic  pressure  associated  with  the  i 

J  mode  of  oscillation,  dimensionless 

q  fluid  velocity  vector 

Re  real  part 

r  plan  polar  coordinate 

r  plane  polar  coordinate,  dimensionless 


Symbol  Description 

r,  horizontal  distance  between  points  (x,y,z) 
and  (£,n,0»  dimensionless 

s  surface  area  of  the  object 

s  surface  area  of  the  object,  dimensionless 

t  time 

T  wave  period 

u  velocity  potential  for  incident  wave,  dimensionless 

u.  velocity  potentials  for  the  radiation  and  scatter 

J  problems,  dimensionless 

X.°  amplitude  of  motion  in  i    direction 


X.°  amplitude  of  linear  or  angular  motion,  dimensionless 

x,y,z  spacial  variables 

x,y,z  spacial  variables,  dimensionless 

C.  i    component  of  force  or  moment  coefficient 

J  associated  with  the  jrn  component  of  oscillation 
of  the  object 

F..  i    component  of  the  dynamic  force  or  moment 

J  due  to  the  j th  mode  of  oscillation 

M. .  added  mass  coefficient 

M.  added  mass  coefficiemt,  dimensionless 

N. .  damping  coefficient 

N  damping  coefficient,  dimensionless 

6  Dirac  delta  function 

£,n,C  rectangular  Cartesian  coordinates,  dimensionless 

9  plane  polar  coordinate 

6.  angular  displacement  of  oscillating  object 

v     =  a  a/g   =   a  tanh  (ah) 
2     2  !/2 

p1      =  (r  +  o 


Symbol  Description 

a  angular  frequency,  a   =    2tt/T 

<p  velocity  potential  for  incident  wave 

<J> .  velocity  potential  for  oscillation  of  object 


J 


or  scatter 


<J> f  total  velocity  potential  for  diffraction  problem 


Subscripts 

i,j         used  to  denote  either  the  mode  or  oscillation 
(or  scatter)  or  the  direction  of  a  force  or 
moment  component,  depending  on  context 
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I.   INTRODUCTION 

The  need  for  the  utilization  of  the  natural  resources 
imbedded  in  coastal  waters  has  led  to  the  need  for  greater 
knowledge  of  the  interaction  of  surface  gravity  waves  with 
large  submerged  objects  such  as  oil  storage  tanks  and  pipe- 
lines.  A  new  problem  arises  in  the  case  of  large  bodies 
immersed  in  shallow  water,  the  scattering  of  the  incident 
wave  and  the  free  surface  effect  become  an  important  considera- 
tion in  the  calculation  of  wave  forces  on  fixed  objects  and 
hydrodynamic  forces  produced  by  oscillation  of  the  body  in 
otherwise  still  water. 

The  analysis  which  considers  these  two  factors  has  come 
to  be  known  as  diffraction  theory.   In  this  approach  separa- 
tion and  viscous  effects  may  be  neglected,  as  was  confirmed 
by  Sarpkaya  and  Garrison  [1],  provided  the  amplitude  of  the 
motion  is  sufficiently  small.   The  problem  is  set  up  in 
terms  of  a  velocity  potential  and  the  potential  function  is 
determined  by  use  of  a  source  distribution.   The  pressure 
on  the  immersed  surface  can  be  found  once  the  source  strengths 
are  known  by  use  of  Bernoulli's  equation  and  the  resulting 
forces  and  moments  are  determined  from  the  pressure  distribu- 
tion by  surface  integration. 

Recently,  Garrison  and  Seetharama  Rao  [2]  worked  out  an 
analysis  for  wave  interaction  with  submerged  objects  based  on 
linear  diffraction  theory  and  applied  the  analysis  to  calculate 
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wave  forces  acting  on  a  submerged  hemisphere.   They  showed 
that  for  small  values  of  the  wave  height  to  sphere  diameter 
ratio,  the  wave  force  coefficients  can  be  well  represented 
by  the  diffraction  theory  and  correlated  as  functions  of  the 
parameters,  2Tra/L  and  h/a  where  h  denotes  the  water  depth,  a 
the  sphere  radius,  and  L  the  wave  length. 

Aside  from  the  wave  force  problem,  an  additional  problem 
may  be  considered  which  is  referred  to  as  the  "radiation 
problem."   This  involves  the  effect  of  a  rigid  body  in  calm 
water  undergoing  forced  harmonic  oscillations  in  its  various 
degrees  of  freedom.   Of  interest  here  are  the  hydrodynamic 
forces  which  may  be  characterized  by  the  coefficients  of 
added  mass,  added  moment  of  inertia,  and  damping.   Although 
often  solved  independently,  the  diffraction  and  radiation 
problems  may  be  solved  concurrently  because  they  are  mathe- 
matically similar. 

Relations  which  eliminate  the  need  to  solve  the  problem 
of  wave  diffraction  past  an  object  by  solving  the  forced- 
oscillation  problem  in  still  water  were  developed  by 
M.  D.  Haskind  in  1957.   They  relate  the  force  and  moment 
components  of  the  diffraction  problem  to  the  asymptotic 
velocity  potentials  for  the  corresponding  radiation  problem. 
That  is,  if  the  waves  produced  by  the  oscillation  of  an 
object  in  still  water  are  known,  Haskind's  relations  may  be 
used  to  obtain  the  wave  forces  (or  moments)  produced  by 
surface  waves  interacting  with  the  same  object  held  fixed. 


12 


Conservation  of  energy  dictates  that  a  balance  must  exist 
between  the  energy  required  to  oscillate  an  object,  and  the 
wave  energy  transmitted  by  surface  waves  across  some  control 
volume  surrounding  the  object  but  at  a  large  radial  distance. 
The  energy  required  to  oscillate  the  object  is  directly 
related  to  the  damping  coefficient,  and  the  damping  coefficient 
may  be  obtained  by  integration  of  the  pressure  over  the  immersed 
surface.   In  view  of  conservation  of  energy  this  same  energy 
appears  at  a  large  distance  from  the  object  in  the  form  of 
radiated  wave  energy.   Therefore,  the  damping  coefficient  may 
be  obtained  by  accounting  for  this  energy  flux.   The  compari- 
son of  this  coefficient  obtained  from  these  two  methods,  i.e., 
the  near  field  and  far  field  solutions,  provides  a  valuable 
check  on  the  accuracy  of  the  numerical  results  generated. 

The  application  of  Haskind's  relations  and  energy  considera- 
tions to  calculate  wave  forces  on  large  submerged  objects  is 
the  purpose  of  this  thesis.   A  computer  program  was  developed 
to  carry  out  the  numerical  calculations  and  to  compare  the 
results  with  those  from  diffraction  and  radiation  theory. 
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II.   THEORETICAL  CONSIDERATIONS 

A.   A  REVIEW  OF  WAVE  FORCES  ON  SUBMERGED  OBJECTS 

Consider  a  rigid  body  of  arbitrary  shape  having  a  charac- 
teristic size  a  submerged  in  water  of  depth  h  as  illustrated 
in  Fig.  1.   It  is  assumed  that  the  amplitude  of  the  motion 
of  the  fluid  is  sufficiently  small  in  relation  to  the  size, 
a,  that  viscous  effects  can  be  neglected.   The  object  may 
or  may  not  be  in  contact  with  the  free  surface  or  the  bottom. 
The  problem  set  forth  herein  deals  with  the  fluid  motion  and 
forces  induced  by  the  small  amplitude  oscillation  of  the 
object  in  its  six  degrees  of  freedom  as  well  as  the  fluid 
motion  associated  with  the  interaction  of  the  fixed  object 
with  a  train  of  regular  waves. 

The  small  amplitude  oscillatory  motion  of  the  body  with 
frequency  about  its  equilibrium  position  is  described  by  the 
relationships 

Xi(t)  =  X°  Re[e"iat]  ,  i  =  1,  2,  3  (la) 

e^t)  =  e!  Re[e'iat]  ,  i  =  4,  5,  6  (lb) 

_  o 

where  X-  denotes  the  amplitude  of  motion  in  translational 

o 

oscillation  and  0.  denotes  the  amplitude  of  motion  in  rota- 
tional oscillation.   The  subscript  i  =  1,  2,  3  corresponds  to 
oscillation  in  the  x,  y,  z  direction,  respectively,  and 
i  =  4,  5,  6  corresponds  to  angular  oscillation  about  axes 
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passing  through  the  origin  of  the  object  parallel  to  the 
x,  y,  z  axes,  respectively.   The  second  problem  considered 
here  is  the  interaction  of  a  train  of  regular  surface  waves 
with  the  object  fixed  in  space.   The  incident  waves  of  wave 
height  H  and  wave  length  L  are  assumed  to  progress  in  the 
positive  x-direction  and  the  fluid  motion  is  assumed  to  be 
incompressible,  irrotational ,  and  harmonic  with  time  dependence 
e     in  all  cases.   Consequently,  a  velocity  potential  exists 
such  that  the  fluid  velocity  vector  may  be  defined  as 

q\  =  Re[V(J)ie"iat]  ,  i  =  1,  2,  ...6  (2) 

where  $ .  denotes  the  velocity  potential  associated  with  the 
motion  induced  by  oscillations  in  the  six  degrees  of  freedom. 
In  the  case  of  regular  wave  interaction  with  the  fixed  object, 
the  velocity  potential  is  given  by  the  sum 

<j>'  =  4>o  +  <f,7  (3) 

where  <J>   denotes  the  velocity  potential  of  the  incident  wave 
in  the  absence  of  the  object  and  <J)7  denotes  the  velocity 
potential  of  the  scattered  wave  due  to  the  presence  of  the 
rigid  body.   For  this  case,  the  fluid  velocity  vector  is 
given  by 

q'  =  Re[V(c^o  +  <$>7)e~10t]       .  (4) 

Having  assumed  irrotationality  of  the  fluid,  <\>  ■    must  satisfy 
the  continuity  equation  which  takes  the  form  of  the  Laplace 
equation,  V2cf>.  =  0  (5) 


16 


within  the  fluid  region. 

Assuming  that  the  amplitude  of  the  motion  is  small,  the 
velocity  squared  terms  in  Bernoulli's  equations  may  be  neg- 
lected.  Therefore,  the  dynamic  pressure  is  given  as 

P.   =  Re[pa4>.e~lcrt],  j  s  1,  2, ...6  (6) 

For  the  second  problem,  involving  wave  interaction  with  the 
fixed  object,  the  pressure  is  given  by 


PT   =  Re[pa(<J>o  +  <J>?)e  iat]      .  (7) 


The  velocity  potential  must  satisfy  certain  boundary 
conditions  in  addition  to  Eq.  (5) .   These  include  the 
linearized  free  surface  boundary  condition, 

-J-  (x,  0,  z)  -  —  <j>.  (x,  0,  z)  =  0  ,  j  =  1,  2,. ..7,         (8) 
9y  g    J 

as  well  as  the  kinematic  boundary  condition  on  the  bottom, 

84). 

—1  (x,  -h,  z)  =  0  ,  i  =  1,  2,. ..7.  (9) 

Since  <J>Q  satisfies  both  Eqs .  (8)  and  (9),  it  follows,  there- 
fore, that  (J)  must  also  satisfy  these  conditions. 

In  addition,  cj>.  must  satisfy  the  boundary  condition  on 
the  surface  of  the  body  defined  by  s(x,  y,  z)  =  0.   This  con- 
dition states  that  for  an  impermeable  boundary,  the  normal 
velocity  of  the  fluid  must  just  equal  the  velocity  of  the 
surface  normal  to  itself.   In  the  case"  of  oscillation  of  the 


rigid  body,  the  kinematic  boundary  conditions  on  the  surface 
of  the  object  take  the  form 


8<t>n 


-  (x,  y,  z,  t)  =  X-n, 


1  x 

9n 


3<J>2  !  o 

-r-  Oc,  y,  z,  t)  =  X?n 
dn  L   y 


3<f>2  *  ° 

-^-  (x,  y,  z,  t)  =  X3nz 
9n 


3*4  •    _    _ 

-T-  (x,  y,  z,  t)  =  0,[(h  +  y)n   -  xn  ] 

dn  4         z     y 


9<J>c 

-3-  (x,  y,  z,  t)  =  95[znx  -  xnz] 


(10) 


^6 
dn 


(x,  y,  z,  t)  =  96[xny  -  (h  +  y)nx] 


where  n  =  in   +  in   +  kn   denotes  the  unit  normal  vector  on 
x   J   y     z 

the  surface  of  the  object  directed  outward  into  the  fluid.   For 
the  case  of  the  interaction  of  regular  waves  with  a  fixed, 
rigid  body  the  normal  velocity  must  be  zero  which  may  be 
written 

d<$>-  d  ({> 

(x,  y,  z,  t)  =  -  — -  (x,  y,  z,  t)  .  (11) 

dn  dn 

Finally,  the  velocity  potential  must  satisfy  the  radiation 

condition  which  permits  only  an  outgoing  progressive  wave  at 

a  large  distance.   Hence  it  is  required  that  <J> .  have  the 
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asymptotic  form  which  satisfies 

-1/2        -  - 
r~   a   -<»    a  fAi"      cosh[k(h+y)l   i(kr,-at)    n     -       r,-. 
c}).(r,  G,  y)  -  A.(0)r1    ■ — - — LJ-J-  e    K      1    ^  -»-  0  as  r,  -►  °°  (12) 

J  3     -1     cosh(kh)  1 

-2-21/2  -- 

j  =  1,  2,  ...7  where  r,  =  (x   +  z  )  '       and  9  =  arctan  (z/x) 

where  A.(0)  is  some  unknown  complex  function  of  8. 

The  wave  number  is  defined  as  k  =  2tt/L  where  L  is  the  wave 

length  and  is  related  to  the  frequency  of  disturbance  according 

to 

|-=  k  tanh  (kh)  .  (13) 

The  velocity  potential  of  the  incident  wave  alone  progres- 
sing in  the  positive  x-direction  is  given  by 

♦  (i,  y,  z)  «  -  ijpl  cosh[k(h_+  y)]  e-ikx         (14) 

cosh(kh) 

_  o 

where  r\      =   H/2  represents  the  amplitude  of  the  incident  wave, 
H  being  the  wave  height. 

In  carrying  out  the  solution  and  to  show  more  clearly  the 
dependence  of  the  solution  on  the  relative  wave  length  parame- 
ter, a  =  27ra/L  =  ka,  and  relative  depth,  h  =  h/a,  a  being  a 
typical  length  of  the  body,  it  is  convenient  to  make  the  space 
variables  and  amplitudes  dimensionless ,  i.e., 

x  =  x/a  ,  y  =  y/a  ,  z  =  z/a  ,  r  =  r/a 

X.  =  X./a  ,  i  =  1,  2,  3   and   X.  =0.  ,i=4,  5,  6 
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and  then  to  introduce  the  pressure  function  u.  by 

o 

ia<J>j(x,    y,    z)/gaX.    =   a   tanh    (ah)    u.    (x,    y ,    z)  ,    j    =    1 ,    2 ,    .  .  .  6 

o 

icr<J>    (x,    y,    z)/gan      =    -au    (x,    y,    z) 

°  (15) 

o 

ia<J>7(x,  y,  z)/gan   =  -au?(x,  y,  z)  . 

The  complex  dimensionless  dynamic  pressure  amplitude  can  now 
be  written  by  use  of  the  linearized  form  of  Bernoulli's 
equation  as 

p.  =  a  tanh  (ah)  u.  (  x,  y,  z)  ,  j  =  1,  2,  ...6 

.    cosh[a(h  +  y)  ]   -iax      r  ■.  ^ 

p'  =  * — * 7-J-L   e     -au7(x,  y,  z) 

cosh (ah) 
where  the  complex  amplitude  of  the  pressure  p  is  defined  as 


3     n  r     -iatn 
«hs-  =  Re[p^  e    ] 

(17) 


PgaX 


P1    n  r  i  -iatn 
— — o  =  Re[p'e    ] 

Pgan 

The  boundary  value  problem  which  describes  the  fluid 
motion  arising  from  the  oscillation  of  the  rigid  body  in  its 
six  degrees  of  freedom  as  well  as  the  scattering  of  the 
incident  wave  can  now  be  written  concisely  in  terms  of  dimen- 
sionless parameters.   The  potential  u.(x,  y,  z) ,  j  =  1,  2,   ...", 
continuous  in  the  fluid  region  is  sought  such  that 

V2u.(x,  y,  z)  =  0   in  y  <  0  (a)  (1?) 
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3u. 

^— L  (x,  0,  Z)  -  a  tanh  (ah)  u.  (x,  0,  z)  =  0      (b) 


8u. 

^— J-  (x,  -h,  z)  =  0   outside  s(x,  y,  z)  =  0       (c) 

dy 


9u. 

3 


(x,  y,  z)  =  h.  (x,  y,  z)   on  s(x,  y,  z)  =  0   (d) 


8n   ^'  ' '      J  j 

,    n    >.    n  rn>.      -1/2  cosh[a(h  +  v)  1   iar     A  ,  , 

u.(r,  6,  y)  -  D.(6)  r   '   \  )  ,  . — =-^-±-  e    i  ->  0  as  r   ->  «   (e) 

j   '   '       1      i      cosh  (ah)  i       v  ' 

where  D.(8)  is  some  dimensionless  complex  function  of  0, 
s(x,  y,  z)  =  0  represents  the  surface  of  the  object  in 
its  undisturbed  position  and  h.  denotes  the  prescribed 
function  which  depends  on  the  mode  of  oscillation  (j  =  1,  2,  ...6); 
the  subscript  7  corresponds  to  scattering  of  the  regular  inci- 
dent wave  due  to  the  presence  of  the  rigid  body.   For  the  six 
degrees  of  freedom  and  scattering,  the  functions  h.  are  given, 
respectively,  by 

h,=n       ,     h-,  =  n       ,     h  _  =  n 
1     x      '      2     y      '      :>     z 

h/i  =  (h  +  y)n   -  zn  ,  hr  =  zn   -  xn  ,  \\,    =  xn   -  (h+y)n      (19) 
4    K     3  J    z  y   5     x     z '   6     y    ^  }  J    x 


h-7  =  w  u\     tn   sinh[a(h+y)l  +  in   cosh[a(h  +  y)l  e 

7    cosh(ah)  ly      i     v      .  ja  x      l  l   ;;J 


lax 


B.   REPRESENTATION  OF  THE  POTENTIAL 

The  solution  to  the  boundary  value  problem  (18)  may  be 
attained  by  use  of  a  Green's  function  having  the  physical 
interpretation  of  a  point  wave  source  of  unit  strength.   These 
sources  are  distributed  over  the  surface  of  the  object 
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according  to  the  source  strength  function  f  so  that  the 
potential  at  some  point  x,  y,  z  within  the  fluid  is  given  by 
the  surface  integral 

Uj  =  **  II     fU>    T]y    C)  G(X'  Yf    ^  *'  n>  °dS       (20) 

s 

where  (£,  n ,  £)  represents  coordinates  on  the  surface  of  the 

object,  f(£,  n,  C)  represents  a  distribution  or  weighting 

-  -2 
function  and  ds  =  ds/a  denotes  the  dimensionless  surface  area 

element. 

The  Green's  function  is  defined,  therefore,  as  the  function 

which  satisfies 

V2G(x,  y,  z;  £,  n,  O  =  6 (x  -  £)  6 (y  -  n)  6(z  -  C)      (21) 

where  6  is  the  Dirac  delta  function,  as  well  as  the  boundary 
conditions  (18b,  c,  e) .   Such  a  function  is  given  by  Wehausen 
and  Laitone  [3]  as 

G(x,  y,  z;  %9    n ,  O  =  \   +  G*(x,  y,  z;  £,  n,  C)      (22) 
where 

1 


G*(x,  y,  z;  £,  n,  C)  = 


R 


+    2p>v#     J        (y   1    v)e"y    tcosh[{i(h   +   nj]    cosh    [y(h   +   y)]  )d 

/  jj   sinn    urn)    -      cosn    (yh)  °o  ^ 

■/   o 

2  2  (23) 

+    i    27T(a      -    v   )    cosh    [a(h   +   n)]    cosh    [a(h   *   y)  ]   j    .. 

a   h   -    v   h   +    v 
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R  =  [(x  -  O2   +  (y  -  n)2  +  (z  -  O2]172 

R'  =  [(x  -  O2   +  (y  +  2h  +  n)2  +  (z  -  c)2]1/2 

r  =  [(x  -  C)2  +  (z  -  O2]172 

2- 

v  =  =  a  tanh  (ah) 

g 

and  P.V.  indicates  the  Cauchy  principle  value  of  the  integral. 
An  alternate  series  form  of  Green's  function  is  also  given 
by  Wehausen  and  Laitone  as 

G(x,y,z;£,n,0  =   2    ~i  cosh[a(h+y)]  cosh  [a  (h  +  n)  ]  • 

a  h-v   h+v 

[YQ(ar)    -    iJQ(ar)]  (24) 

+    4   U       -2 2 cos[yk(h+y)]cos[yk(h+n)]Ko(ykr) 

k=l  y,h-v  h+v 

where  J   and  Y   denote,  respectively,  Bessel  functions  of  the 
o       o        >     r        /  » 

first  and  second  kind  of  order  zero  and  K   denotes  the  modi- 

o 

fied  Bessel  function  of  the  second  kind  of  order  zero.   The 
quantities  y,  are  the  real  positive  roots  of  the  equation 

yk  tan  (yRh)  +  v  =  0  .  (25) 

The  solution  to  the  boundary  value  problem  now  rests  on 
the  determination  of  the  source  strength  function  f.   Taking 
the  normal  derivative  of  the  potential  u.  and  applying  boundary 
condition  (18d)  yields  the  following  integral  equation  from 
which  f  is  to  be  determined: 
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4tt  JJ     fj(^,n,C)|§(x,y,z;?,n,Ods  =  h.(x,y,z)  ,  j  =  1,  2, ...7  (26) 


where  3G/8n  may  be  obtained  by  straight  forward  differentiation 
of  (22)  and  (24). 

C.   NUMERICAL  SOLUTION 

A  numerical  procedure  can  be  devised  by  approximating  the 

actual  contour  by  a  profile  of  a  finite  number  of  facets. 

Proceeding  in  this  direction,  the  surface  of  the  object  is 

partitioned  into  N  area  elements  of  size  As.  where  the  sub- 
1  3 

script  takes  on  values  j  =  1,  2,  ...N.   Since  the  source 
strength  function  f,  occurring  in  the  integrals  in  (20)  and 
(25),  is  a  continuous,  well-behaved  function,  these  integrals 
may  be  approximated  by  the  following  summations: 

f  .  a.  .    =  h  .  (27) 

nj   ij     ni  v      J 

u  .  =  f  . B .  .  (28) 

ni     ni  ij  v  ' 

where 


1    ff        9G  ,  ., 

a..  =  -. —  //    tt—  (x  .  ,y  •  ,  z  -  ;x  .  ,y  .  ,  z  . )  d: 
13    4tt  JJ         8n  ^  l'7!'  1'  j9/j'    3 


(29) 


As. 
3 


6-.  =  i-  //    G(x. ,v. ,z. ;x. ,v. ,z.)ds.  (30) 

ij    4tt  JJ  K    i'-i'  1'  j  3       3 

As 


1 


th 


In  these  expressions,  u  .  denotes  the  potential  at  the  i 
nodal  point  on  the  object  associated  with  either  the  n  ,  mode 
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of  o   illatio/i,  or  scattering  (n  =  7)  of  the  incident  wave. 

The    L«g  rations  are  carried  out  over  the  finite  surface 

are-  -?leir.     As.  and  f.   denotes  the  value  of  f  at  the 

cent  :oil  points  of  The  area  element.   The  approximation  made 

heTc  ds  cons:  Lent  in  that  as  N  approaches  infinity,  the 

appTcximatio: v,  given  in  (27)  and  (28)  approach  the  exact 

forms . 

■••■'-■.'"-nee     :.-a\         es  a,  3,  and  h  are  numerically  evaluated, 

(2'  yean     SLy  b(  inverted  by  use  of  a  standard  computer 

subroutine  t   'btain  the  solution  for  f.   u   is  then  obtained 

n 

from  (29).   There  are,  however,  certain  difficulties  connected 
with  the  eva" action  of  (29)  and  (30)  which  must  first  be  con- 
sidered.  Specifically  for  the  special  case  when  i  =  j,  the 
1/R  and  o(l/K)/8n  terms  are  singular  as  R  approaches  0,  and 
secondly  in  "be  .valuation  of  (29)  and  (30)  due  to  the  singu- 
lar natuic  ojl  [l/(y  tanh  (uh)  -  a  tanh  (ah)]  at  y  =  a  occurring 
in  the  infinite  integral  in  (23)  and  its  normal  derivative. 
These  difficulties  are,  however,  circumvented  by  subtracting 
■  ou'i  :the  sing  laxity  and  carrying  out  its  integration  analytical' 
ly.   Thes   d  fyils  are  carried  out  in  some,  yet  unpublished, 
notes  c.t  -("..  '■',    Garrison. 

D.   HYDRODYNAMIC  FORCES  AND  MOMENTS 

The  forces  and  moments  caused  by  the  dynamic  fluid  pressure 
acting  upon  the  immersed  surface  of  the  object  may  be  obtained 
fror,  the  integrals, 
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F^Ct)  =  -  (1  or  £)  fTpjhidi   ,   i,  j  =  1,  2,  ...6    (31) 

Fi(t)  =  -  (1  or  a)  JJP,hids    ,   i  =  1,  2,  ...6       (32) 

where  F.  denotes  the  i    component  of  wave  force  (or  moments) 

and  F..  denotes  the  i    component  of  forces  arising  from  the 

j    component  of  body  motion.   The  coefficient  1.0  is  used  for 

the  case  of  a  force  (i  =  1,  2,  3)  whereas  a  is  used  when  F 

denotes  a  moment  (i  =  4,  5,  6). 

For  purposes  of  presentation  of  the  numerical  results, 

dimensionless  force  coefficients  are  defined  as 

F.       .  x 
r  imax    io.  -,   0   ,  /,,•> 

Ci  =  — -_  e   l   ,   1  =  1,2,3  (33) 

Pga  n 

and 

F  i  =  1,  2,  3 

C.  =  _iiH|2L_  ei6ij  =  -M..  -  iN..,     .    -   ?     ,  ^34) 

in      -4V°      J  in     ii'     i  =  1,  2,  ...6. 

J         pga  X.  J  J  J  y       > 

The  corresponding  expressions  for  the  moment  coefficient  are 

n  imax    io.  .   r  r  r-zc^ 

Ci  =  — ^T^"  e   x    ,i  =  4,5,6  (35) 

Pga  n 

and 

F.       .-  i  =  4,  5,  6 

=  _2JM*_  el5ij  =  -  M..  -  iN-.  ,  (36) 

1J    ^!  1J      13     j  -  1,  2,  ...6. 

■f-Vi 

The  complex  coefficients  C.  relate  to  the  i  component  of 
wave  force  (or  moment)  coefficient  while  C-.  denotes  the  i 
component  of  force  (or  moment)  coefficient  associated  with 
the  j    component  of  oscillation  of  the  object.   The  phase 
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shift  angles  6.  and  5..  relate  the  phase  of  the  force  to  the 
crest  of  the  incident  wave  and  displacement  of  the  body, 
respectively.   The  real  and  imaginary  parts  of  the  dimension- 
less  force  coefficients,  M. .  and  N..,  are  called  the  added 
mass  and  the  damping  coefficients,  respectively. 

Using  (31)  and  (32)  in  conjunction  with  (16)  and  (17) 
and  the  definitions  (33-36),  the  force  coefficients  may  be 
written  in  dimensionless  form  as 

C.  =  jj[auy(x,y,z)  -  C0^(^]]]  h.(x,y,z)ds ,   i=  1,2, ...6  (37) 
s 

and 

Cij  =  J I  Uj (x>y>z)hi(x>y>z)ds   »   i,  j  =  1,  2,  ...6.       (38) 


Once  u.(x,y,z)  is  obtained  from  (28),  the  coefficients  can  be 
obtained  from  (37)  and  (38)  by  numerical  quadrature. 


E.   HASKIND'S  RELATIONS  AND  ENERGY  CHECK 

Although  it  may  be  presupposed  that  the  numerical  solution 
outlined  here  will  converge  upon  increasing  the  number  of 
partitions,  modern  computers  are  still  limited  by  storage 
capacity  and  computer  time  is  expensive.   In  view  of  this, 
it  is  necessary  to  keep  the  partition  size  as  large,  i.e., 
the  number  of  partitions  as  small,  as  accuracy  considerations 

permit.   The  determination  of  the  effect  of  the  partition 
size  on  accuracy  becomes,  therefore,  an  important  considera- 
tion in  order  that  practical  limits  may  be  established.   One 


27 


available  method  is  the  comparison  of  the  numerical  results 
with  analytical  results  where  closed  form  solutions  exist. 
This  approach,  however,  is  limited  to  a  few  simple  shapes; 
for  more  general  configurations  no  such  check,  of  course, 
exists . 

A  method  which  is  valid  for  checking  the  soundness  of 
the  numerical  results  of  any  shape  is  the  use  of  the  so- 
called  Haskind's  relations  as  well  as  an  energy  balance,  both 
of  which  are  developed  in  the  Appendix.   The  former  is  used 
to  verify  the  solution  for  the  diffraction  problem  and 
the  latter  the  solution  of  the  radiation  problem.   Using  the 
asymptotic  form  of  Green's  function  as  given  in  (24)  in  con- 
junction with  (20)  the  following  relationship  for  the  damping 
coefficient  is  so  obtained 


n..  -  i-  [  a 

11    2p  L 


v 


? 


a "  h  -  v 


f77  ff 


ur  ru  ,  \i  - iapn cos  (3-9 i  j 
:osh[a  (h  +  n)  ]e    ^1    ^    d: 


'de 


(39) 


where  p,=  (£2  +  C2)     and  $  =  arctan  U/O  .   This  relates 
the  damping  coefficient  to  the  far  field  behavior  of  the  solu- 
tion.  It  is  in  effect,  an  energy  balance  between  the  energy 
required  to  oscillate  the  body  and  the  energy  transmitted 
across  some  control  surface  a  large  distance  from  the  body. 

A  relationship  somewhat  similar  to  (59),  known  as  Haskind's 
relations,  may  be  obtained  for  the  wave  force  (or  moment) 
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coefficient  and  relates  the  wave  produced  at  infinity  by  the 
body  oscillating  in  the  j    mode  to  the  j    component  of 
wave  force  as 


Cj    cosh  (ah)  J  J     f  ,(£,n,0 


s   J 


cos 


ur    r-u,    m  -  iap-,  cos  (3-tt)  j  (40) 

h[a(h+n)]e    Kl    *■    yds.  v  J 


This  result  represents  the  Haskind's  relation  for  finite  depth 

as  given  by  Seetharama  Roa  [4] . 

Equations  (39)  and  (40)  represent  relations  for  the 

damping  and  wave  force  (or  moment)  coefficient  based  on  the 

behavior  of  the  far  field  solution.   A  comparison  of  these 

results  with  N.   and  C.  obtained  from  an  integration  of  the 
13       3  5 

pressure  over  the  immersed  surface,  i.e.,  as  obtained  from 
(37)  and  (38)  ,  provides  a  convenient  and  valuable  self-check 
on  the  accuracy  of  the  numerical  results.   These  results  are 
not  limited  to  special  configurations  and  may  be  applied  to 
arbitrary  shapes.   Equation  (40)  is,  however,  limited  to 
symmetry  with  respect  to  the  x-y  plane,  a  condition  which  is 
generally  satisfied  by  practical  shapes. 
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III.   COMPUTER  PROGRAM 

In  order  to  feasibly  utilize  the  method  described  in 
section  D  of  THEORETICAL  CONSIDERATIONS,  a  computer  program 
was  developed  to  carry  out  the  indicated  calculations.   Although 
accuracy  was  a  controlling  factor,  as  it  should  be,  the  pri- 
mary requirement  was  to  hold  the  computer  run  time  to  a 
reasonable  value.   To  achieve  this,  every  advantage  was  taken 
to  employ  the  symmetry  in  the  equations  and  to  eliminate  redun- 
dant computations  by  generating  certain  constants  and  matrices. 

The  program  consists  of  several  subroutines,  each  solving 
a  specific  part  of  the  problem.   This  allowed  the  development 
of  the  program  to  be  simplified  and  the  use  of  it  was  made 
easier  in  that  only  a  small  main  program  was  required  to 
execute  the  program.   Subroutine  GEODAT  reads  the  input  geo- 
metrical data,  generates  the  matrix  h  as  defined  in  (19)  ,  and, 
in  conjunction  with  GCOEFF,  calculates  certain  geometrical 
parameters.   GARRIS,  the  primary  subroutine,  then  generates 
the  a  and  3  matrices  by  calling  either  of  two  subroutines. 
The  first  (GREEN)  calculates  G  and  3G/8n  based  on  the  integral 
form  given  in  (22)  and  (23)  whereas  the  second  (GREENS)  makes 
evaluations  on  the  basis  of  the  alternate  series  formulations 
given  in  (24).   For  elements  of  the  a  and  8  matrices  corres- 
ponding to  small  values  of  (ar)  GREEN  is  utilized,  while  for 
large  values  of  (ar)  GREENS  is  used.   With  the  exception  of 
the  diagonal  elements,  the  majority  of  the  elements  of  a  and 
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Figure  2:   Program  Flow  Chart 
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$  are  calculated  by  the  latter  subroutine.   This  situation 
proved  to  be  most  fortunate  because  the  series  form  converges 
rather  rapidly  and,  consequently,  requires  much  less  computer 
time  than  the  integral  form. 

In  order  to  check  the  validity  of  the  two  subroutines, 
values  of  G  were  computed.   A  comparison  between  GREEN  and 
GREENS,  Figs.  3,  4,  5}  for  constant  a  and  decreasing  r,  shows 
that  at  small  values  of  a,  the  two  methods  give  similar 
results  whereas  at  intermediate  values  of  a,  the  GREEN  solu- 
tion becomes  very  erratic  at  larger  values  of  r.   At  high 
values  of  a,  however,  the  two  subroutines  render  similar 
solutions  despite  the  higher  frequency  of  oscillation.   This 
stems  from  the  fact  that  the  amplitude  of  oscillation  decays 
extremely  fast. 

Subroutine  COMAT  then  determines  the  source  strength 
function  f  by  solving  the  complex  matrix  equation  (27). 
Once  the  matrix  f  has  been  generated,  the  potential  function 
u  is  calculated  by  (28).   The  pressure  at  points  on  the 
immersed  surface  as  well  as  the  resulting  forces  and  moments 
are  then  determined  in  subroutine  FORCES. 

Finally,  the  Haskind's  relations  and  the  energy  check  are 
employed  to  calculate  the  wave  force  and  moment  coefficients 
and  the  damping  coefficients  in  COCHCK.   Here  also,  a  compari 
son  is  made  between  these  values  and  those  found  in  FORCES. 
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Figure  3:   Comparison  of  Equations  (22)  and  (24);  a  =  0.2 
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Figure  4:   Comparison  of  Equations  (22)  and  (24);  a  =  1.0 
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IV.   DISCUSSION  OF  RESULTS 

A.   VERTICAL  CYLINDER 

In  order  to  study  the  accuracy  of  the  computer  program, 
numerical  results  were  generated  for  both  the  wave  force 
problem  and  the  problem  involved  in  the  oscillation  of  the 
object  in  otherwise  still  water.   Simple  geometric  shapes, 
including  vertical  circular  cylinders  and  submerged  spheres, 
were  studied  and  compared  with  closed  form  solutions  where 
they  exist. 

As  indicated  in  Section  II,  the  accuracy  of  the  numeri- 
cal results  is  dependent  on  the  size  of  the  grid.   Garrison 
[Ref.  5],  for  one,  found  that  the  fineness  of  the  grid  size, 
N,  must  be  increased  as  a  increases.   This  results  from  the 
fact  that  the  kernels  of  the  integral  equations  oscillate 
rapidly  as  the  parameter  a  increases. 

A  series  of  computer  runs  was  made  for  a  vertical  circu- 
lar cylinder  resting  on  the  bottom  and  piercing  the  free 
surface.   Two  runs  consisting  of  cylinders  with  grid  sizes 
N  =  120  and  N  =  252  were  made  for  a  relative  depth  of  h  =  2.0 
in  order  to  determine  the  effect  of  the  grid  size.   Another 
was  made  with  N  =  120  and  h  =  6.0  to  determine  the  effect 
of  depth.   A  fourth  run  was  made  with  N  =  240  and  h  =  4.0, 
but  this  time  the  grid  was  divided  into  t     :zes.   The 
bottom  half  of  the  cylinder  had  N  =  96  with  an  element 
aspect  ratio  (width/height)  of  0.525  whereas  the  top  half  had 
N  =  144  and  an  aspect  ratio  of  0.785.. 
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Figure  6:   Definition  Sketch 
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Within  the  range  of  practicality  of  a,  a  =  0  to  a  =  3.0, 
the  numerical  results  for  the  horizontal  wave  force  coeffi- 
cient and  the  coefficient  of  moment  about  a  bottom  line  axis 
compared  favorably  with  the  closed  form  solution  given  by 
MacCamy  and  Fuchs  [Ref  6],  Figs.  7  and  8.   These  results 
show  that  within  this  range  of  the  parameter  a,  there  is  such 
a  small  deviation  in  the  values  obtained  for  the  different 
grid  sizes  that  the  two  curves  are  superimposed.   As  the 
frequency  is  increased,  however,  it  was  found  that  deviation 
from  the  closed  form  results  increased  more  rapidly  for 
the  coarser  grid  than  the  finer  one.   It  is  also  note- 
worthy that  the  magnitude  of  the  moment  coefficient  relative 
to  the  wave  force  coefficient  is  small  at  low  frequencies 
and  increases  as  the  frequency  increases.   This  result  is 
as  expected  since  at  small  values  of  wavelength,  wave  action 
is  concentrated  near  the  free  surface  whereas  at  large 
values  of  wavelength,  wave  action  is  felt  at  deeper  depths, 
thus  decreasing  the  moment  arm. 

Figures  9  and  10  represent  the  hydrodynamic  force  and 
moment  coefficients  associated  with  the  cylinder  oscillating 
in  still  water.   These  show  that  for  zero  frequency,  the 
force  coefficient,  C,,,  approaches  approximately  2tt  and  4tt 
for  the  h  =  2.0  and  h  =  4.0  cylinder,  respectively.   At 
h  =  6.0  a  value  of  approximately  6tt  was  obtained.   When  this 
is  converted  to  an  added  mass  coefficient  the  result  is  1.0. 
This  is  as  expected  since  at  zero  frequency  the  free  surface 
acts  as  a  rigid  boundary  and  the  flow  corresponds  to 
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Figure  9:   HYDRODYNAMIC  FORCE  COEFFICIENT  FOR 
A  VERTICAL  CIRCULAR  CYLINDER 
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two-dimensional  flow  past  a  circular  cylinder.   The  classical 
value  for  the  added  mass  coefficient  in  this  case  is  1.0. 

B.   SPHEROID 

A  series  of  runs  for  varying  a  was  made  for  a  half 
submerged  sphere.   One  consisted  of  a  grid  size  of  200 
points  with  the  relative  depth  h  =  5.0  and  the  others  were 
made  with  a  grid  size  of  264  points  and  h  =  5.0  and 
h  =  10.0.   Havelock  [Ref.  7]  determined  the  added  mass  and 
damping  coefficients  for  a  floating,  heaving  sphere.   His 
results  were  later  confirmed  by  Kim  [Ref.  8]  who  also 
obtained  such  values  for  floating,  surging  (or  swaying) 
spheroids.   Although  the  results  found  by  Havelock  and  Kim 
were  for  infinite  depth,  the  values  found  in  this  study  for 
h  =  5.0  and  10.0,  Figs.  11-16,  compare  quite  favorably. 
There  is  apparently  little  depth  effect  beyond  h  =  5.0 
except  at  small  values  of  a  where  it  is  noted  the  curves 
corresponding  to  h  =  5.0  and  h  =  10.0  differ. 

Kim  [Ref.  9]  also  determined  the  horizontal  and  vertical 
wave  force  coefficients  as  did  Milgram  and  Halkyard  [Ref.  10]. 
Comparison  of  these  values  with  those  calculated  in  this 
study,  Figs.  15  and  16,  show  good  agreement  with  the  former, 
but  some  difference  was  noted  with  the  latter,  especially 
for  the  vertical  force  coefficient  where  the  value  of 
Milgram  and  Halkyard  was  almost  three  times  greater  at  a  =  2.0 
It  would  appear  the  Milgram  and  Halkyard  results  are  in  error. 
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For  a  fully  submerged  spheroid  there  were  only  two 
available  checks  known.   One  is  that  the  value  of  the  added 
mass  coefficient  for  a  sphere  in  an  infinite  fluid  is  0.5. 
Figure  18  shows  that  for  the  numerical  procedure  used  in 
this  study  a  value  of  0.51  was  obtained  for  a  relative  depth 
of  5.0.   The  second  test  involves  the  use  of  Haskind's  rela- 
tions which  relates  the  vertical  force  on  an  object  fixed 
in  a  monochromatic  wave  train  with  the  heave  damping 
coefficient  N_~  such  that 

m*     a     sinh  (2ah)     r    2  rA1, 

NL  =  T  ^-r- — 7 ■  -I       [TV  C~>  (41  ! 

22    2  2 ah  +  sinh  I  2 ah)   2  v  J 

where  N*.,  denotes  the  damping  coefficient  in  heave.   A 
comparison  of  the  damping  coefficient  found  in  this  study 
with  that  using  (41)  is  shown  in  Table  1.   Using  the  results 
of  the  diffraction  theory  as  the  standard,  the  percentage 
deviations  from  it  for  the  Haskind's  relations  are  also 
shown.   The  comparison  shows  that  for  values  of  a  greater 
than  0.5,  all  of  the  results  obtained  were  within  4%  of  the 
taken  standard.   For  values  of  a  less  than  0.5,  the  differ- 
ences were  considerable  and  this  was  believed  to  be  due  to 
the  low  value  of  the  relative  depth  h.   A  similar  compari- 
son was  made  for  the  264  point  floating  sphere  using  h  =  10.0 
The  results  are  shown  in  Tai  "   2.   Since  the  present  results 
satisfy  Haskind's  relations,      ay  be  concluded  the  Milgrai 
and  Halkyard's  results  are  de  finitely  in  error. 
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Table  1:   Damping  Coefficient  in  Heave  for  a  Submerged  Sphere 

N  =  200        h  =  5.0 


a 

N22 

N* 
1N22 

%      deviation 

0.20 

0.0267 

0.0199 

25.46 

0.30 

0.1033 

0.0914 

11.52 

0.50 

0.3382 

0.3252 

3.84 

0.80 

0.4433 

0.4329 

2.34 

1.00 

0.3745 

0.3698 

1.25 

1.50 

0.1610 

0.1589 

1.30 

2.00 

0.0514 

0.0515 

-0.19 

3.00 

0.0033 

0.0034 

-3.03 

4.00 

0.0002 

0.0002 

0.0 

5.00 

0.0 

0.0 

0.0 

Table  2:   Damping  Coefficient  in  Heave  for  a  Floating  Sphere 

N  =  264        h  =  10.0 


a 

N22 

N* 
1N2  2 

%   deviation 

0.20 

0.5650 

0.5278 

6.58 

0.40 

0.7071 

0.6919 

2.15 

0.60 

0.6968 

0.6920 

0.69 

0.80 

0.6250 

0.6071 

2.86 

1.00 

0.5356 

0.5202 

2.87 

1.20 

0.4525 

0.4388 

3.02 

1.40 

0.3799 

0.3679 

3.15 

1.60 

0.3176 

0.3075 

3.18 

1.80 

0.2654 

0.2555 

3.73 

2.00 

0.2211 

0.2304 

-4.20 
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C.   HASKIND'S  RELATIONS  AND  ENERGY  CHECK 

As  previously  indicated,  important  checks  on  the  numerical 
technique  and  the  final  results  can  be  made  by  application  of 
Haskind's  relations  and  the  energy  check  which  are  developed 
in  the  Appendix.   These  checks  were  applied  to  all  computer 
runs  made  and  the  percent  difference,  using  the  surface 
integration  of  the  pressure  results  as  the  standard,  is 
plotted  in  Figs.  19  and  20. 

In  general,  it  was  noted  that  dependence  of  the  results 
on  the  grid  size,  N,  is  much  greater  for  the  forces  obtained 
by  a  surface  integration  of  the  pressure  when  calculated  by 
use  of  Haskind's  relations.   It  can  be  seen  from  Figs.  19 
and  20  that  as  frequency  decreases,  the  percent  deviation 
decreases.   At  large  values  of  a  the  error  becomes  quite 
large.   The  effect  of  decreasing  the  grid  size  or  increasing 
the  depth  is  also  clearly  shown.   For  the  120  point  grid  the 
errors  are  much  larger  than  for  the  240  point  grid  in  the 
case  of  the  h  =  2.0  circular  cylinder.   Similarly  the  errors 
for  the  120  point  grid  are  much  smaller  for  h  =  6.0  than 
for  h  =  2.0. 
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V.   CONCLUSIONS 

A  computer  program  based  on  linear  wave  theory  has  been 
developed  to  calculate  wave  forces  acting  on  submerged 
objects  of  arbitrary  shape.   In  order  to  reduce  the  computa- 
tion time,  the  restriction  of  symmetry  with  respect  to  two 
planes  passing  through  the  center  of  the  object  was 
incorporated. 

On  the  basis  of  the  results  presented  in  Section  IV, 
the  following  conclusions  appear  justified: 

1.  All  of  the  checks  and  comparisons  were  successful 
and  therefore  it  appears  that  the  present  method  yeilds 
accurate  results  up  to  approximately  a  =  3.0  depending  on 
the  grid  size. 

2.  The  horizontal  force  and  moment  coefficients  in  all 
cases  increase  with  the  wavelength  parameter  a,  at  first, 
reach  a  peak  and  then  decrease  as  a  increases. 

3.  The  vertical  force  coefficient  for  a  floating 
spheroid  starts  with  a  value  of  it,  corresponding  to  zero 
size  to  wavelength  ratio,  and  decreases  rapidly  at  first 
and  later  slowly  as  a  increases. 

4.  Haskind's  relations  and  the  energy  check  give 
valuable  checks  on  the  accuracy  and  show  the  effect  of 
increasing  the  grid  size  or  the  depth. 
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APPENDIX 
HASKIND'S  RELATIONS  AND  ENERGY  CHECK 

A.   HASKIND'S  RELATIONS 

The  dimensionless  wave  force  and  moment  components,  F-(t), 
on  an  object  of  arbitrary  shape  may  be  obtained  by  substituting 
for  P'  from  (16)  and  (17)  in  (32)  and  simplifying.   For  the 
diffraction  problem,  this  gives 

re  9u.        . 

F.(t)  =  Re[a  JJ   (u?  +  uq)  ^  ds  e~iat],  i  =  1,  2,  ...6   (A.l) 

s 

where  u.  denotes  the  radiation  potentials. 

l  v 

The  scatter  potential  can  be  eliminated  from  (A.l)  through 

the  application  of  Green's  reciprocal  theorem,  applied  to  G 

and  u.  as 
3 

\\\     [u.(x,y,z)V  G(x,y,z;£,n,C)-G(x,y,z;£,n,C)V  u.(x,y,z)] 
R 


dxdydz  =   /  /  [u.  (x,y ,z)-^-(x, 

Sr+S,+S   +S 

f   b   °° 


lj^»/»*J^n 


8u. 

- G  (x  ,y  ,  z  ;  %  , n  ,  ?)^_i(x  ,y ,  z)  ]  ds  , 

in  the  region  R,  as  shown  in  Fig.  21.   Since  an  incompressible 
fluid  was  assumed,  both  u7  and  u.  satisfy  the  Laplacian  in 
in  the  region  R,  Eq .  (18a), 
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ff  3u7     3u. 

J  /  (u.~— -  -u7T-i)ds  =  0 

J  J_ .         v  i3n    73n  ^ 


(A. 2) 


s  +  s  +s,  +s 
f  b 


The  contributions  to  the  surface  integral  in  (A. 2)  from  the 

surfaces  Sr,    S,  and  S   vanish  since  u-  and  u-  satisfy  the 

f»   b       °°  7       l       J 

free  surface  and  bottom  boundary  conditions  and  the  radiation 
condition  at  R   =  °°.   Also,  from  (18d)  in  conjunction  with  (14) 
and  (15),  (A. 2)  yields 

ff  3u.  ff  3u.  3u 

/  /    (u7+u    )    ~  ds    =     /  /    (u  »~i   -u.    ,j-£)    ds    .  (A. 3) 

JJ    v    7      o7    3n  JJ     v    o3n  l    3n    J  K 

s  s 

Again  using  Green's  theorem,  the  surface  integral  on  the  right 
hand  side  of  (A. 3)  may  be  replaced,  thus  giving 

r  r  3u.         r  r  3u.      3u 

/  /  (U  +u  )    *  ds  =  -  /  /   (u   ,-i  -u.  3-^)  ds  .       (A. 4) 
//    7   o^  3n         J  /   v  o  3n     l  3n  J  y 


s 

00 


By  using  the  radiation  condition  (18e)  and  since  x  =  r  cos 
it  can  be  shown  that 

r  r  3u.      3u  r  r 

j   /   (u  —i  -u.  -s-2)  ds  =  -   J  /   u  u.[iar(l-cos  6) 
JJ      v  o  3n    i  3n  J  ^    J  J  o  iL    v 


s  s 

oo  oo 


(A. 5) 


\~]    ds. 
2rJ 


Substituting  for  u   and  u.  from  (14)  and  (15)  and  then 
&       o       1 


simplifying 
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3u.  3u 


if      (u      Tnr-   "u-    -3 )    ds    =  - * 

//       v    o    °n  1    <m    '  ,2,   ,  ^ 

J  J  a   cosh    (ah) 


rsinh(2ah)        h,      iar 
[~ TiE "  +   2  J    e 

2tt 


J  Di(9)r'1/2[l/2-iar(l    -    cos    6)]    eiar   cos    G    de  (A. 6) 


-|      O    -y*  ^    y-i    ^  Q 

Inasmuch  as  r  -*■  °°  and  e  is  a  fluctuating  quantity, 

the  real  part  of  (A. 6)  equals  zero.   Thus 

2tt 

/t.  ra^-       -l/2>i        Q^   iar  cos  0  ,a 
D.(6)iar   '  (1  -  cos  6)  e  d6 

o 

27T 

.   -1/2   iar  /    D.(9)(l  -  cos  9)  e-iar(1  "  cos  9)d9   (A'7) 
=  -iar   '   e    J  r  M         ' 


Since  ar  >>  1,  the  solution  to  the  integral  on  the  right 
hand  side  of  (A. 7)  may  be  obtained  by  the  method  of  stationary 
phase  (Stoker  [11]).   Hence 


/ 


n  r«wi        n-\  -iarfl  -  COS  9)  Jc. 

D. (9)  (1  -  cos  9)e     v         '    d0 


o 

ar  ■> 


(A. 8) 


0   1/2  .  ,ir   0   >. 
=  2C— )    el(T  "2ar)  D.(tt) 


and  after  simplifying 
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// 


3u.'         3u  ,      1/2  , 


odn        1    dn  ^a    7  ,2,   ,  . 

cosh    (ah) 

00 

[£inhi2ahl+  eiJD.(lF)  (A. 9) 

z  a  l 

From    (A.l),     (A. 3),    (A. 4)    and    (A. 9),    simplification   gives 

1/2 
F.  (t)    =    Re{-ii2jLi) fSinMJahl  +   h]    ei  Cj-ot)^  (ff)  }  (AaQ) 

1  coshz(ah)     za  x 

and  taking  the  modulus  on  both  sides  gives 

1/2 
.  il^a) [Sinh(2ah)  +  h]  ,  M  ,  >  .  .  ^  2j_fi     (An) 

cosh  (ah) 

In  actuality,  (A. 11)  represents  six  different  equations 
which  are  known  as  "Haskind's  relations."   They  relate  the 
i    component  wave  force  (or  moment)  coefficient  to  the  waves 
produced  at  infinity  by  the  object  oscillating  in  the  i   mode. 

B.   ENERGY  CHECK 

+■  v> 
Equation  (31)  represents  the  i    component  of  the  dynamic 

force  (or  moment)  arising  due  to  the  j    mode  of  oscillation 

of  the  object.   In  the  case  of  the  radiation  problem,  they 

may  each  by  expressed  as  the  sum  of  two  components,  one  in 

phase  with  the  acceleration  and  the  other  in  phase  with  the 

velocity  of  the  body  in  the  form 

F.  . (t)  =  -  M..aX.(t)  -  N..aX.(t)  ,  i,j  =  1,  2, ...6         (A. 12) 

where  M. .  represents  the  added  mass  (or  added  moment  of 
inertia)  and  is  defined,  when  compared  to  (34),  (36)  and  (37), 
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as 

-3     C  T  i  =  l,2,3 

M^  =  -pa   Re[//  u^(x,y,z)hi(x,y,z)ds]  >    *    =i      2,       .6    (A-13a) 

s 
for  a  force  and  for  a  moment  as 

-4       r  r  1  =  4,5,6 

Mi;j  =  -pa^  Re[JJ  u.(x,y,z)hi(x,y,z)ds],   ..  =  ^  2#>_#6j  (A. 13b) 

s 
and  N. .  denotes  the  damping  coefficients  (in  relation  to 

energy  transport)  and  is  defined,  when  compared  to  (34),  (36) 

and  (37)  ,  as 

i  =  1,  2,  3 


paa   Im[  f  J  u^  (x,y,z)hi(x,y,z)ds]  ,   j  =  i   2    .6   (A. 14a) 


N.  . 

"  s 


for  a  force  and 

i  =  4,  5,  6 


N 


-4        r  r  1=4,^,0 

i.  =  -paa   Im[J  J  u . (x,y , zjhu  (x,y , z)ds] ,   j  =  j   2    .6   (A-14b) 


for  a  moment . 

The  dimensional  force  (or  moment)  component,  in  the  i 

direction,  on  the  oscillating  body  due  to  oscillations  in  the 

i   mode  may  then  be  obtained  from  (A. 12)  by  equating  j  to  i. 

Since  the  forces  and  moments  exerted  by  the  body  on  the  fluid 

are  equal  and  opposite  to  these,  the  energy  transmitted  by 

the  object  to  the  surrounding  fluid  during  one  period  is  given 

by 

T 

E.  =  -f  (1  or  a)  F..  (t)  L(t)  dt  (A. 15) 

where  a  is  used  for  a  force  (i  =  1,  2,  3)  and  the  factor  1  for 
a  moment  (i  =  4,  5,  6).   After  carrying  out  the  integration, 
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(A. 15)  may  be  rewritten  as 


E..  =  7rpazabX?zNi:L,   i  =  1,  2, ...6  (A. 16) 

where  N. .  denotes  the  dimensionless  damping  coefficient  which 
is  defined  for  a  force  as 

N. .        i  =  1,  2,  3 

v  '  To?    •   J  ■  i.  2>---6 

and  for  a  moment  as 

N. 


N 


ij        i  =  4,  5,  6 

ii      -4 
J    paa       j  =  1 ,  2 , . . . 6 . 


If  n°  is  the  amplitude  of  the  outgoing  progressive  waves 
at  S   due  to  oscillation  of  the  body  in  the  i   mode,  then 

00  '  ' 

this  amplitude  can  be  related  to  the  asymptotic  velocity 
potential,  u. ,  by  using  the  dynamic  free  surface  boundary 
condition  and 


n°(r,9)  =  aX?  tanh(kh)  a  |u.(r,8,0)|  .  (A. 17) 

The  total  energy  transmitted  over  one  period  can  now  be 
written 

-4  4 
~    Lim  pL  r,       2kh    -,  a  a   v02  n 
h.  =      ^—  [1  +  J  a.   r 

1    f->«>         sinh(2kh)    g     * 


I 


2tt 

|u. (r,e,0) |2d6  ,  i  =  1,  2, ...6.  (A. 18) 

o 


Now,  by  comparing  (A. 16)  and  (A. 18),  the  damping  coeffi 
cients  may  be  written  in  terms  of  the  asymptotic  velocity 
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potential ,    u . ,  .as 


Lim    1    ...__,_  ,._i_A-ri  2ah 


N..    =    ":J"  4  tanh(ah)[l    +      ■    ZTi  u] 
2  v.      ^  L  smn  (2 an) J 


/ 


2tt 


ui(r,6,0)  |2d6  ,   i  =  1,  2, ,..6.  (A. 19) 


C.   NUMERICAL  EVALUATION  OF  HASKIND'S  RELATIONS  AND  ENERGY  CHECK 

The  right  hand  sides  of  (A. 11)  and  (A. 19),  involving  the 
asymptotic  velocity  potentials,  are  solved  by  the  use  of  the 
distribution  functions  f.,  obtained  from  the  numerical   solu- 
tion  of  the  radiation  problem,  in  conjunction  with  the  series 
form  of  Green's  function  (24). 

For  r,-*-00  and  y  =  0,  considerable  simplification  is  possible 
in  the  asymptotic  form  of  G  and  may  be  written  as 

G(x,  0,  z  ;  £,n  ,  O  =  -  — j± ~ cosh(ah)  cosh  [a  (n  +  h)  ] 

a  h-v  h+v 


2_   1/2  ei(ar1-iT/4) 


(— — )     exvail  "'^  .  (A. 20) 


Th 


us  the  asymptotic  potential,  u.,  is  given  by 

u.(r,  e,  0)  =  -  4(-?-)1A  (^-a2)cosh(ah) 
iv  '   *  J  2Ki\aJ  2,   2, 

a  h-v  h+v 


// 


1/2    .  ,         7T. 

£•(£,  H,  O  cosh  [a(n  +  h)]f1    ellarl"  TJds  .  (A. 21) 


s 

Since  it  is  large,  r,  may  be  replaced  by  r  and  moved  outside 
the  integral.   Further,  as  shown  in  Fig.  22,  for  large  r 

r,  =  r-p,  cos  (8-6) 
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2  2    1/2 

where    p,    =    (£      +    C    )         ,    and  u.    may  be   rewritten   as 


Ui(r,     6,  0)    =    -i(2    a) 


-1/2    .2      2.         ,  e   ,.       -1/2    .,        tta 
(v    -a    )cosh(ah)  /       i(ar-T) 

27      2T~ 
a   h-v   h+v 


JJ   f±U,    n,  c)cosh[a(n+h)]  eiaplcos  (6_e)  ds  . 
s 

Substituting  (A. 22)  into  (A. 19)  and  simplifying  gives  the 

result 


(A. 22) 


,      2   2 
N,,  =  *  [   a  "v 


ii    Ztt  l  2,       2\ 

a  h-v  h+v   o 


v   °      l/2s 


C)cosh[a(n+h) ] 


.e-iap1cos(B-6)  +  e-iap1cos (6+0)^d, 


de  . 


(A. 23) 


Once  u.(r,  6,  0)  is  known,  D-(tt)  may  be  obtained  by  using 
the  radiation  condition  (18e) .   Thus 

1/2 


lar 


DiO)  =  ui(r,TT,0)r    e"~x  (A. 24) 

and  after  simplifying  by  means  of  trigonometric  identifies, 
(A. 11)  becomes 

Ci  =  cosh  (ah)  J  J         fi(?,v,C)cosh[a(n  +  h)]e"ia^lcos(6"^ds,  (A. 25) 

l/2s 

where  i  =  1,  2,**6  since  for  an  axisymmetric  body,  the  only 
cases  of  interest  are  those  pertaining  to  i  =  1,  2,  and  6. 

A  comparison  of  the  results  obtained  from  (A. 23)  and  (A. 25) 
with  those  calculated  by  integration  of  the  pressure  over  the 
submerged  surface,  Eqs .  (37)  and  (38),  are  a  valuable  self- 
check  on  the  accuracy  of  the  numerical  results. 
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